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Abstract. 

Approximate analytical solutions in closed form are obtained for the 5-dimensional 
Bohr Hamiltonian with the Woods-Saxon potential, taking advantage of the Pekeris 
approximation and the exactly soluble one-dimensional extended Woods-Saxon 
potential with a dip near its surface. Comparison to the data for several 7-unstable 
and prolate deformed nuclei indicates that the potential can describe well the ground 
state and 71 bands of many prolate deformed nuclei corresponding to large enough 
“well size” and diffuseness, while it fails in describing the /3i bands, due to its lack of 
a hard core, as well as in describing 7-unstable nuclei, because of the small “well size” 
and diffuseness they exhibit. 
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1. Introduction 

The advent of critical point symmetries BE!, related to shape/phase transitions in 
nuclear structure M, has stirred interest in analytical solutions (exact or approximate) 
|5j of the Bohr Hamiltonian [ 6 j. In addition to the infinite well potential, used in 
the critical point symmetries E(5) [Tj and X(5) [ 2 ], related to the transition from 
spherical to 7 -unstable (soft with respect to triaxiality) [TJ and to prolate deformed 
jE] nuclei respectively, solutions involving the Davidson j9[ IT). HI, ECS :l3j, Kratzer 
[m da usi 1123 , and Morse unman]] potentials (shown in figure 1 ) have been given. 
When applied to the bulk of nuclei for which energy spectra and B(E2) transition rates 
are experimentally known [ 21 ] , all these potentials provide good and quite similar results 
[T 2 l T3l flly jl 6 l ITT1 I9l 120] . despite having quite different shapes. 

The question arises if the success of the above mentioned potentials is due to the 
form of the Bohr Hamiltonian alone, or if there are potentials which, when plugged into 
the Bohr Hamiltonian, will not be able to provide satisfactory results for nuclear data. 

I 11 this work we consider the Woods-Saxon potential [22], which has been 
extensively used in nuclear physics in a different context, namely as a single-particle 
potential [23]. The potential, shown in figure 1(c), reads 

u ws(P) = - e2a( ^° o) + 1 , (1) 

where Uo, a and /?o are non-negative free parameters. This study is motivated by several 
reasons. 

1) For a —* 00 the WS potential reduces to a finite square well potential, which 
has also been used in relation to critical point symmetries [23, [25], in addition to the 
infinite well potential BB- 

2) I 11 contrast to the Davidson [9], Kratzer [IT] , and Morse [TS] potentials, which 
possess a hard core, the WS potential does not have a hard core. We will find out the 
consequences of this difference. 

3) The Woods-Saxon potential is also interesting from the mathematical point of 
view, since it has no closed analytical solution for the spectrum, even in one dimension 
and for vanishing angular momentum [261127] . Only an analytical wavefunction is known 
in one dimension and for vanishing angular momentum [ 27] . 

First we manage to produce approximate solutions in closed form for the WS 
potential with a centrifugal barrier, exploiting its resemblance (after applying the Pekeris 
approximation [28] to it) to a modified spherically symmetric WS potential [29] known 
to possess exact analytical solutions [26]. Taking advantage of these solutions within the 
framework of the Bohr Hamiltonian we subsequently fit several 7 -unstablc and prolate 
deformed nuclei, pointing out the successes, but also the failures of the WS potential. 

Numerical calculations involving generalized forms of the Bohr Hamiltonian have a 
long history. 

1) Extensive early numerical calculations for vibrating axially symmetric deformed 
nuclei have been performed in the framework of the Rotation Vibration Model [30l 0T] . 
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Figure 1 . (Color online) The Kratzer [T4] (a), Davidson j9j (b), Woods-Saxon [22] 
(c) and Morse R8] (d) potentials, for special values of their free parameters, used for 
simplicity. All quantities shown are dimensionless. 


2) Numerical solutions have subsequently been provided using a general form for 
the potential energy [32, [33], giving emphasis on the restrictions imposed on the form 
of the potentials by symmetry constraints. 

3) Extensive numerical results using a general collective model employing general 
forms of both the kinetic energy and the potential energy have been obtained initially 
[ 34 ] by using a basis provided by Hecht [35], and subsequently [36] 137] by exploiting the 
basis of a 5-dimensional (5D) harmonic oscillator [38] [39]. 

4) The recent clarification [40] of the group theoretical structure of the Davidson 
potential, when used in the framework of the Bohr Hamiltonian, led to the development 
of the Algebraic Collective Model [CL, 02] 03] 03], which allows the efficient numerical 
calculation of spectra and transition probabilities of nuclei of any shape. 

In addition to the y-unstable and the prolate deformed nuclei, the description of 
triaxial nuclei within the framework of the Bohr Hamiltonian was very early attempted 
[45, 46] and is still attracting considerable attention [25, @7] 0E, 49] [50]. 

On the other hand, studies on the microscopic foundation of the Bohr Hamiltonian 
have a long history, some early efforts reported in [51] [52] [53]. Efforts towards the 
derivation of the Bohr collective Hamiltonian through the adiabatic approximation of 
the time-dependent Hartree-Fock-Bogolyubov theory (the ATDHFB method), as well 
as through the generator coordinate method with the Gaussian overlap appoximation 
(the GCM+GOA method) have been reviewed in [54]. Recently, the derivation of the 
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parameters to be used in a 5D Hamiltonian for quadrupolc vibrational and rotational 
degrees of freedom from Relativistic Mean Field (RMF) calculations [55] [56, 5?J [58j is 
acquiring momentum. 


2. Modified Woods—Saxon potential with a centrifugal barrier 


2.1. Pekeris approximation for the centrifugal barrier 


The Pekeris approximation [28] is a well known method for finding approximate 
analytical solutions for potentials involving exponentials. It has been introduced 
in relation to the Morse potential [T8] with a centrifugal barrier in studies of 
rotational-vibrational spectra of diatomic molecules [28]. The basic idea of the Pekeris 
approximation is to rewrite approximately the centrifugal term using exponentials 
resembling the ones appearing in the rest of the potential, with the aim to be able 
to “absorb” the centrifugal term into the potential. 

This approximation has been used in relation with the Morse potential in diatomic 
molecules [28. .59], as well as with the pseudo-centrifugal term for Dirac particles within 
Morse potentials [60]. Furthermore, a supersymmetric improvement of the Pekeris 
approximation in the case of the Morse potential has been worked out ED- Numerical 
results reported in these papers indicate that the approximation works well, especially 
for relatively low-lying vibrational and rotational states (see in particular the detailed 
results reported in [6T] ). 

The Pekeris approximation has further been used in relation to various exponential- 
type potentials [[62], to the Wei Hua oscillator [63], to the Rosen-Morse and Manning- 
Rosen potentials [63], as well as with the pseudo-centrifugal term for Dirac particles 
within Rosen-Morse potentials [65]. In these cases, various functions appear in the 
original potentials, as well as in the Pekeris approximation terms. 

The spherically symmetric Woods-Saxon potential in three-dimensions with a 
centrifugal barrier reads 

m = £W0 )+Km = - e ^-w + 1 + (2) 

where Uws stands for the Woods-Saxon potential, U c stands for the centrifugal term, l 
is the angular momentum, while Uq, a and /?o are non-negative free parameters. 

It is known [26] that this potential cannot be solved exactly. However, closed 
analytical solutions can be derived using the Pekeris approximation [28]. As it was 
mentioned above, the basic idea of the Pekeris approximation is to rewrite approximately 
the centrifugal term using exponentials resembling the ones appearing in the rest of the 
potential, with the aim to be able to “absorb” the centrifugal term into the potential. 

In the present case we wish to write the centrifugal term in the approximate form 


UM = 


i(i + i) 


Cn + 


c 1 


+ 


c 2 


e 2 a (/ 3 - p 0 ) + 1 ' [ e 2a(/3—/3 0 ) + ]_]2 


/3 2 




(3) 
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In order to use a more compact notation, we introduce the change of coordinates 
x = (j3 — Po)/Po- Then the exact centrifugal term takes the form 

5 


UJx) = 


where 


5 = 


(1 + x) 2 ’ 

1(1 + 1 ) 

Pi ' 


Expanding this function into binomial series one gets 
U c (x) = 5(1 — 2x + 3x 2 — Ax 3 + ...)• 

On the other hand, using the same notation, the approximate form of 

Ci , C 2 


U c (x) 


Co + 


+ 


e 2a/3 0 x _|_ |g2 a/3 0 x _|_ 1 J2 

Expanding the exponentials in Taylor series we get 

Ci , C 2 \ Po 

Xo 


U c (x) 


c ° + Y + t)~ it (Cl + ° 2)x + mf 2 * 1 + 


where 


o o — 


2 a' 


Equating the coefficients of equal powers of x in dHJ) and (J 8 ]) we get 


n , 4a 0 12 al 

c °- l -J^ + YP 


n _ 8a 0 

C ‘~H 


48cio 

’ 


Using the approximate expression of d3J) in 


U(P) » 5C 0 - 


(C 0 - 5C,) 


+ 


Po 

3), we obtain 

6C 2 


c = 48 «o 


PI 


(4) 

(5) 

( 6 ) 

becomes 

(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 


g2a(/3—/3 0 ) 1 ' [ e 2a(/3-/3 0 ) + 1]2’ 

where indeed the centrifugal term has been “absorbed” in the potential. 

The accuracy of approximating the potential of (J 2 ]) by (fill) will be discussed in 
detail in section 5, in relation to the parameter values occurring for the nuclei under 
consideration (see figs. 4, 5, 6 and relevant discussion in section 5). 


2.2. Modified spherically symmetric Woods-Saxon potential 

On the other hand, the modified spherically symmetric Woods-Saxon potential, which 
presents a dip near its surface |29j, has the form 

U 0 We 2a C~M 


V(P) = 


( 12 ) 


e 2a(/3-/3 0 ) _|_ X [ e 2a(/3-/3 0 ) + 1]2 ’ 

where Vo, W, a and Pq are non-negative free parameters. This potential, shown in figure 
2 , is known [26] to possess exact solutions of the one-dimensional Schrodinger equation 
in terms of Jacobi polynomials [ 66 ] 


+ n (P) oc 


1 


e 2a(/3-/3 0 ) _|_ X 


6/2 


1 


g2a(^-^ 0 ) q. X 


c/2 


p(6,c) 


(13) 
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Figure 2. (Color online) The Woods-Saxon potential with a dip near its surface, 
given in is shown for special values of its free parameters, used for simplicity. All 
quantities shown are dimensionless. 


where 


with 



c 


W 

Pn = \ 1 H-o ~~ n + 1)’ 


and the condition [ B 6 | 



b,c> — 1 


has to be fulfilled, implying the restriction 
b + c > - 2 . 





The strengths of the potential of (fT21) can be expressed in the form [26] 
V 0 = a 2 (b 2 — c 2 ) > 0, 

W = a 2 [{b + c)(b + c + 4n + 2) + 4 n{n + 1)] > 0, 


(14) 

(15) 

(16) 

(17) 

(18) 
(19) 


which imply stronger restrictions on b, c. Indeed, from (TT9ll for n — 0 one should have 
{b + c)(b + c + 2 ) > 0 , which implies either b + c < — 2 , which is not allowed by (fT7|h or 

b + c> 0 , ( 20 ) 


which is a stronger restriction. Higher values of n do not impose any further restrictions, 
since (fT9l) leads to b + c < — 2 (n + 1 ), which is not allowed by ffl7l) . or b + c > — 2 n, 
which is weaker than (JTTJ). It should be remembered that the quantity W, given in 
(fT9]h cannot vanish, since in that case the potential of (fl2|l would collapse to the usual 
Woods-Saxon potential, which is known to have no closed analytical solutions for its 
spectrum [26, 27]. 

The eigenvalues of the energy are known to be 

x V*\ Vo 
2 ' 


E n~ 4 


Pn + 


« 4 Pn 


( 21 ) 
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Using the identity 

g2a(/3—/3 q) 


|g2a(/3—/3 0 ) _|_ 1]2 e 2a(/3-/3 0 ) + 1 [g2a(/3-/3 0 ) _|_ l]2 ’ 

the potential of (1T2li can be rewritten as 

Vq + W W 


V(P) = 


+ 


e 2a(P-p 0 ) + 1 [ e 2a(/3-/3 0 ) + 1 ] 2 ‘ 

2.3. Spectrum of the modified Woods-Saxon potential with centrifugal barrier 
We remark that this potential has the same form as the potential of (TTT]h with 


( 22 ) 


(23) 


leading to 


V 0 + W = U 0 - SC U 

* = «,- 4 ' (; + 1) 


W = 5C 2 , 

121(1 + 1 ) 


W = 


aPo " a 2 P £ 

Substituting these expressions in (121 p and taking into account the first term in dTIT) 
the eigenvalues of the energy for the potential of (fill) are found to be 


(24) 


(25) 


E n = 


1(1 + 1 ) 
PI 


12 a| 

PI 


Pn 

4a 0 


U 0 a 0 


8i(Z+l)ag 

4o 3 


Pr 


Uo 

2 ’ 


where 


192Z(Z + l)an , 

Pn = \0 + - „ 4 - ( 2 « + !)• 


(26) 


(27) 


Therefore we have managed to obtain the energy eigenvalues of the WS potential 
with a repulsive barrier, as approximated through the Pekeris approximation, by 
exploiting the known solutions of the modified spherically symmetric WS potential, 
without having to solve the Schrodinger equation anew. 

We remark that the condition of ()20l) leads to severe restrictions for n. Indeed, 
using m we have 

b + c = p n > 0, (28) 

which, using (1271) . leads to 


1 

n < - 
~ 2 


-l + \/l + 


192Z(Z + 1)oq 


Pi 


(29) 


Thus in the case of / = 0, only n — 0 is allowed. 

This result is valid in a 3-dimensional space. Following the same path one can see 
m I68J that the result is also valid in D dimensions, with the angular momentum l 
replaced by 

l D = l + -7 ^. (30) 

In the special case of a 5-dimensional space, l is replaced by 

L = Z + 1. 


(31) 
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3. Bohr Hamiltonian for 7 -unstable nuclei 


3.1. Energy eigenvalues 
The original Bohr Hamiltonian [6] is 
h 2 r 


// - — 


2 B 


+ - 1 - 

/^d/r d/3 (3 2 sin 37 (Ty 


d d 

— sm 37— 

07 


4/1 2 


£ 

fc=1,2,3 


Ql 


(7 - Ink ) 


Sill" 


+ m 7 ), 


( 32 ) 


where f3 and 7 are the usual collective coordinates describing the shape of the nuclear 
surface, Qk {k = 1, 2, 3) are the components of angular momentum, and B is the mass 
parameter. In what follows we are going to consider h 2 = 2 B = 1. 

Assuming that the potential depends only on the variable /3, i.e. V(/3, 7 ) = U(/3), 
one can proceed to separation of variables in the standard way U EL using the 
wavefunction 


* 03 , 7.00 = 


(33) 


where 6, : (i = 1,2,3) are the Euler angles describing the orientation of the deformed 
nucleus in space. 

I 11 the equation involving the angles, the eigenvalues of the second order Casimir 
operator of SO(5) occur, having the form A = t(t + 3), where r = 0, 1 , 2 , ... is the 
quantum number characterizing the irreducible representations (irreps) of SO(5), called 
the “seniority” [69]. This equation has been solved by Bes P3). 

The “radial” equation yields 


d 2 f 


U(P) + 


(r+ l)(r+ 2 ) 


HP) = E n , T f(/3). 


d(3 2 V P 2 

The effective potential appearing in this equation coincides with that of 
formal replacement of l by r + 1. Substituting in 


E = 


where 


(t + 1)(tT2) 


PI 


1 + 


12ar 


Pi 


Pv 


4 a n 


U 0 ao 


(34) 

), with the 
we obtain the energy eigenvalues 
,a\ 2 


8(r+l)(r+2)a5 


Pn 


Pn,r ~ \ / 1 + 


192(r + l)(r + 2)ag 


P 4 o 


- (2n+ 1). 


-T- (35) 


(36) 


The restriction of ([29]) in this case reads 
-1 + 


1 

n < - 
~ 2 


192(t + 1)(t + 2)oq 


1 + 


K 


(37) 


As we shall see in subsection 5.2, only n = 0 is acceptable for the parameter values 
occurring in real nuclei. The violation of the condition of (l37]i even for n — 1 has as a 
consequence to push the n — 1 bands too low in energy, making them unphysical. 


















Analytical solutions for the Bohr Hamiltonian with the Woods-Saxon potential 


3.2. Rescaling the energy eigenvalues 
Eq. (l25l) takes the form 

8a 0 (r + 1 )(t + 2) 


U 0 = V 0 + 


PS 


( 38 ) 


implying that if Vo is a constant, then f / 0 depends on r, and vice versa. We are going 
to consider Vo as being the constant quantity. Then we can use this equation in order 
to eliminate Uq from (1351) . Using the notation 

A = ao/ Po, 

the equation for the energy becomes 

En,r = (T + + ^ (1 + 12 » ~4i) - - igg ~ 

Po 16 «6 Pn,r 2 

with 


(39) 


(40) 


i>n,T= y 1 + 192(r + 1)(t H- 2)A i — (2n + 1). (41) 

In order to reduce the number of free parameters, we fit nuclear spectra leaving out 
overall scales, as it is done in the E(5) and X(5) models [HE], Instead of fitting the raw 
experimental levels El, where L is the angular momentum, the quantities 

(49) 

£/ 2 — H/q 

are used. In other words, the energy of the ground state is subtracted from all levels, 
and then each level is divided by the energy of the first excited state, which in even-even 
nuclei is the energy of the first excited state with L — 2. The connection between L 
and r is described in detail in subsection 13.41 

In the present case this has the following consequences on (140]) : 

1) The last term, —Vo/2, plays no role, since it is a constant and cancels out. 

2) Using the rescaling 


Po — 


Po 


a o — 


ao 


(43) 


%/vT VvT 

a common factor Vo appears in all terms of (1401) . which therefore cancels out when 
above is used. 

3) Using A = a 0 /p 0 = a 0 /p 0 , a common factor (3% can be taken away from the 
denominator of all terms of ( 1401) . 

Then the equation to be used for the energy fits becomes 


E ntT — (t + 1)(t + 2)(1 + 124 2 — 44) — 


Pn,T 

164 2 


A 2 P A 0 

pI,t 


(44) 


where p n , T is still given by (1411) . 

This is the equation used in the fits, from which the parameters A and Po can be 
determined. Then a 0 is calculated from a 0 = Ap 0 . 
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The Woods-Saxon potential is essentially different from zero within its radius, which in 
this case is (3q. 

One way to fix the scale is to identify the radius of the potential, /3 0 , with the (3 exp 
value obtained from the experimental values of B(E2] Of —> 2f) [Tlj- Then from (1431) 
one simply has 


K,= 


Pi 



(45) 


Notice that in this case ao can be calculated directly from a o = Af3 exp . 

It should be recalled at this point that while the B(E2]0f — > 2^) values are 
experimental quantities which are model independent, the deformation parameter is 
model dependent mi. In the present case, a uniform charge distribution out to the 
distance R(6, <p) is assumed, with zero charge beyond it, while the nuclear radius is 
taken to be i?o — 1-2 A 1//3 fm. Then the deformation is connected to the transition 
probability by 


P = 


An 

3 ZRl 


B(E2;0f -> 2{ 


(46) 


where B(E2] 0j 




2+) is measured in e 2 b 2 . It should be noticed that f3 is 


the deformation parameter defined in the Bohr framework [ 6 j. A slightly different 
deformation parameter e is used in the framework of the Nilsson model [72], connected 
to the Bohr parameter by e = 0.946/3 (see p. 125 of [72], or eq. (2.82) of [73]). 


3.4- The spectrum 


The spectrum is characterized by the 0(5 )dS 0(3) symmetry [69j[70j. r and L are the 
quantum numbers characterizing the irreps of 0(5) and SO(3) respectively. The values 
of angular momentum L contained in each irrep of 0(5) (i.e. for each value of r) are 
given by the algorithm HU 

T = 3z/ a + A, Z/ A = 0,1,..., (47) 

L = A,A+1,...,2A —2,2A (48) 


(with 2A — 1 missing), where zz A is the missing quantum number in the reduction 0(5) 
D SO(3), and are listed in Table 1. 

The ground state band (gsb) has n — 0 and levels L g — 0, 2, 4, 6, ..., for which 
zz A = 0 and L g = 2r. Thus within the gsb one has 

(r+l)(r + 2)= (L + 2) 4 (L + 4) . (49) 

The quasi -71 band has n — 0 and levels with L 7 = 2, 3, 4, 5, ..., which are 
characterized by zz A = 0 and L 1 = 2r — 2 for L 1 being even, or by L 1 = 2r — 3 for L 1 
being odd. Therefore they exhibit the following degeneracies 


2=4 


3 7 — 4 7 — 5 7 — 6 7 — 


’9i 


7 7 = 8 7 = 10„,... 


(50) 
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The quasi-/3i band has levels Lp — 0, 2, 4, 6, .... There are two choices: 

1) n = 0, in which the levels have ua — 1 and Lp = 2r — 6, leading to the following 
degeneracies 

Op = 6 g, 2 p = 8g, 4 p = 10 3 , 6 p = 12g, . . . (51) 

2) n — 1, in which the levels have va = 0 and Lp = 2r. Then no degeneracies of 
this kind occur. 


4. Bohr Hamiltonian for prolate deformed nuclei 


4-1. Energy eigenvalues 


If the potential has a minimum around 7 = 0, as it is the case for prolate deformed 
nuclei, the angular momentum term in (l32]i can be written [2] as 


3 


sin 2 

k= 1 


i^ k 27ri \ ~ o (Ql + Ql + Q 3 ) + Q 
(7 - ~k) 3 


sin 2 7 


(52) 


Exact separation of /? from the rest of the variables can be obtained for potentials of 
the form Bi 


U(P, 1 )=u(P) + AA. (53) 

Concerning the 7 and Euler angles, the solution provided by Bes cm is not valid any 
more, thus we seek, as is customary in the literature (see Chapter IV of the review 
article [5] for further details) wave functions of the form 

= hW < 54 ) 


where 0j (j = 1, 2, 3) are the Euler angles, T>(9j ) represents Wigner functions of these 
angles, L stands for the eigenvalues of the angular momentum, while M and K are the 
eigenvalues of the projections of the angular momentum on the laboratory-fixed 2 -axis 
and the body-fixed V-axis respectively. In these wave functions the dependence on the 
Euler angles is entering through the Wigner functions, as in the many examples reviewed 
in [5], while K is assumed to be a good quantum number, an assumption which is not 
a priori justified, since strong K mixing can be present. 

The 7 -equation occurring from this separation of variables has been solved in [12], 
using a potential 


U7) = (3c)V, 


(55) 


its eigenfunctions written in terms of Laguerre polynomials and its energy eigenvalues 
given by 


e 7 — (3 C)(n 7 + 1), 
while the separation constant is 

x _ 


C — 2c, n 7 = 0,1,2,3, 


(56) 


(57) 
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It should be noticed at this point that the 7 potential used in the Bohr Hamiltonian has 
to be periodic in 7 j 6 l |33l [75]. The potential of (l55]h used in several solutions 
is a lowest order approximation, representing the first 7 -dependent term in the Taylor 
expansion of the proper periodic potential. There has been recently extensive work 
involving proper periodic 7 potentials in the Bohr Hamiltonian [5] 08] [76j ES 78( [79, 


The radial equation takes the form 

1 d o4 d L(L + 1) A 
p±dp^ dp + 3/3 2 + / 32 +m ^) 

Transforming into Xl by the relation 

AT (A) 


UP) — ^l(P)- 


we obtain 


UP) = 


x"M + 


P 


2 > 


L(L+ 1 ) 


3 +X+ —~u(P) 


P 2 


Xl{P) = 0 . 


(58) 


(59) 


(60) 


We remark that this equation looks very similar to (1341) . with L ^ L , +1 ' > + A + 2 replacing 


(r + l)(r + 2). Therefore from (l35l) we see that the energy eigenvalues become 


E n ,L = 
where 


('zyn+q + A + 2 j 


A 


1 + 


12 ai 


A, 


o f Ua 8(^i)+A + 2)ag\ 
Pn.L \ I U ° a ° W 0 


4a n 


Pn,L 




Pn,L — \ / 1 + 192 


L(L + 1) 


2 + A + 2 ) -77 — (2 n + 1). 


Ao 4 




(62) 


The restriction of in this case reads 


1 

n < - 
~ 2 


L(i + 1) + A + 2 )| 


-I + 4 /I+ 192 


(63) 


As we shall see in subsection 5.1, only n = 0 is acceptable for the parameter values 
occurring in real nuclei. The violation of the condition of (l63|l even for n = 1 has as a 
consequence the failure of the relevant fitting attempts. 


4-2. Rescaling the energy eigenvalues 

Following the same rescaling procedure as in the 7 -unstable case, the equation used in 
the fits becomes 


^-(^ +A+2 ) (1+12i2 - 4i) "S «u’ 




where 


Pn.L — 1 /1 + 192 I — 1 ———! + A + 2 ) A 4 — (271 + l). 


(64) 


( 65 ) 
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This is the equation used in the fits, from which the parameters A, /? 0 , and C 
(appearing in A) can be determined. 


V 



Figure 3. (Color online) Effective potentials for L = 10 for some Yb isotopes, obtained 
from (l23l) . The parameters are taken from Tabled while W is given by (l25l) with Z(^+l) 
replaced by L H+ 1 ) + \ _|_ 2, where A = 3C from (l57l) and (l56l) . The quantities shown 
are dimensionless. 


5. Numerical results 


5.1. Prolate deformed nuclei 


As a first test, the spectra of 63 nuclei with R 4 / 2 > 2.9 have been fitted. The quality 
measure 


a N 


Ei=i ( E i( ex P) - Ei(th)f 
(lV-l)(£(2+))2 


( 66 ) 


has been used in the rms fits. 

The results for 46 nuclei are shown in Table [21 The following comments apply. 

1) Only the ground state band and the quasi -71 band have been included in the 

fits. 

2) When trying to include the quasi-/3i band in the fits, trying to correspond it to 
the n — 1 case, the fits fail. The failure is due to violation of the condition of ([63]) . 
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The parameters obtained from fitting the experimental spectra of 162 ~ 176 Yb, shown 
in Table [21 exhibit the following features. 

1) The depth of the potential, Vo, exhibits a clear minimum near the middle of the 
neutron shell (l^Ybiozi), where maximum deformation is observed, as indicated by the 
clear maximum exhibited by /3q = fd exp , the experimental values of deformation obtained 
from the B(E2; 0 + —» 2 + ) [71] . 

2) The highest diffuseness do occurs nearest to the shell closure, while the lowest 
appears near the middle of the shell. 

3) In retrospect, the selection of the parameters used in the fits was physically 
meaningful, fio is related to the “well size” d. while A is related to the ratio of the 
diffuseness parameter a 0 over the average width of the potential /3 0 . 

4) The potentials obtained from the Yb isotopes are shown in figure 3, corroborating 
the above made remarks about maximum depth and diffuseness nearest to the shell 
closure, and minimum depth and diffuseness near the middle of the shell. We remark 
that the dip near the surface of the modified spherically symmetric WS potential is 
becoming very large, dominating the overall shape of the potential. It is instructive to 
compare the shapes of the potentials appearing in figure 3 to these in figure 2, since 
they come from ([23]) and (fl2l) respectively, which are equivalent. 



Figure 4. (Color online) Exact (dashed lines) and approximate (solid lines) effective 
potentials for L = 0, 6, 12, 18, 24, 30 for 232 Th, obtained from (J2J and (1111) respectively. 
The parameters are taken from Table O providing the quantities needed from (f5[) . ((Oil . 
and m while Uq is given by (l25l) with l (l + 1) replaced by L H+ 1 ) -p A -f- 2, where 
A = 3C from (1571) and (15^1) . The corresponding energy levels, determined from (1711 . 
are also shown. The quantities shown are dimensionless. See subsection 15. II for further 
discussion. 

The spectra obtained for some rare earths and actinides are shown in Tabled Good 
agreement with the data is obtained. 

A question arising at this point is the accuracy of the Pekeris approximation used 
in the derivation of the relevant formulae. In order to check this, we show in figure 4 
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both the original exact potential (dashed lines), obtained from ([2]), and the approximate 
potential (solid lines) occurring after the Pekeris approximation, obtained from (fill) , for 
six different values of the angular momentum L , in the case of 232 Th. In addition, the 
relevant energy levels, obtained from (j6Tj) , are shown. The following comments apply. 

1) Near their bottoms, the exact and the approximate potentials nearly coincide. 
They exhibit minima very close to each other, while in addition the two potentials form 
similar wells around the minima. 

2) At higher energies, the approximate potential wells are slightly displaced to the 
left, in relation to the exact potentials. However, at the energy of main interest, i.e. at 
the energy of the relevant level, the widths of the two potentials are nearly equal. 

3) The approximation appears to be better at higher values of L, i.e. at higher 
energy levels. This is expected, since the fits were obtained through an rms procedure. 


u u u v 



Figure 5. (Color online) Same as figure 4, but for L = 0, 6, 12, 18, for 170 Yb (a) and 
176 Yb (b). See subsection 15.II for further discussion. 

The same observations hold for all nuclei shown in Table [2j with two exceptions, 
1(2 Yb and 176 Yb, the latter exhibited in figure 5. While 170 Yb (and the rest of the Yb 
isotopes shown in Table [2]) exhibit the qualitative behaviour described above, in 176 Yb 
(as well as in 7 ' 2 Yb) we remark that the approximation breaks down, since for the 
lowest values of L the relevant energy level lies higher than the corresponding barrier 
on the rhs of the exact potential. In other words, we still get fits of good quality, which 
are obtained using the approximate potential, but this potential is not similar to the 
exact potential any more, preventing us to draw any conclusions regarding the exact 
potential. It should be noticed at this point that the approximate potentials (solid 
lines) for Yb isotopes shown in figures 3 and 5 are identical, up to a displacement by 
a constant. Indeed, the potentials of figure 3 come from (1231) . while the potentials of 
figure 5 come from (fill) , the latter being displaced in relation to the former by the term 
5Cq. This term is included in figure 5, since it is necessary for the comparison to the 
exact potentials (dashed lines) given by (|2]) . 

1,2 Yb and 176 Yb are not isolated cases. Similar breakdowns of the approximation 
have been found in 17 more nuclei, shown in Table [4j In all cases, good fits of quality 
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comparable to those shown in Tableware obtained, but they correspond to approximate 
potentials being very different from the exact ones. 

From the numerical results obtained for all the nuclei mentioned above, we see that 
the breakdown of the approximation occurs for “well size” /3 0 < 1-8 and/or diffuseness 
do < 0.43. Since /3 q appears in the denominator of the rotational term, this result is 
understandable. The lower /3 0 becomes, the more important the rotational term, in 
which the Pekeris approximation has been used, becomes, thus the approximation to 
the exact potential deteriorates. 

We therefore conclude that good fits are obtained for all nuclei studied, but only if 
the “well size” /3o fulfils the condition j3o > 1.9 and the diffuseness fulfils the condition 
do > 0.44 the approximate potentials are quite similar to the corresponding exact 
potentials. 



Figure 6. (Color online) Exact (dashed lines) and approximate (solid lines) effective 
potentials for r = 0, 3, 6 (L = 0, 6, 12) for 118 Xe, obtained from <[2j) and (fill) 
respectively. The parameters are taken from Table [5j providing the quantities needed 
from ©, ©, and m, while W is given by (l25l) with/(Z + l) replaced by (t + 1)(t + 2). 
The corresponding energy levels, determined from (1401) . are also shown. The quantities 
shown are dimensionless. See subsection 15.21 for further discussion. 


5.2. 7 -unstable nuclei 

The spectra of 38 nuclei with i ? 4 / 2 < 2.6 [ 98 > 10 o>io 2 r Uj io 2 ,io 4 ,io 6 ,io 8 ,no,ii 2 ,ii 4 ,ii 6 p c |^ 
108,110,112(7^ 118,120,122, 124 ,126,128, 130 , 132 , 134 -^-g 130 , 132 , 134 p a 134 , 138 Q e 142 Qd 156 Er 

186 , 188 ,190,192,194, 196 , 198 ,200p(.j j iave been fitted, using the same quality measure. 

The following comments apply. 

1) Only the ground state band and the quasi -71 band have been included in the 

fits. 

2) When trying to include the quasi-/3i band in the fits, trying to correspond it to 
the n = 1 case, one finds it lying too low in energy, due to the violation of the restriction 

of prp. 

3) When trying to include the quasi-/?! band in the fits, trying to correspond it 
to the n — 0 case, the levels of the /h-band in most cases are overestimated. There is 
nothing to be done against this, since the degeneracies of (j5TT) fix the position of the 
theoretical levels. Thus this failure is due to the fact that the experimental data do not 
exhibit exactly the 0(5) degeneracies. 
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The parameters obtained from fitting the experimental spectra of 118 - 134 Xe, shown 
in Table [5j We remark that in all cases the “well size” /3 0 is much lower than 1.8, while 
in addition the diffuseness is much lower than 0.43 . The same holds for all of the 
38 nuclei considered. Therefore the approximate potentials bear little relevance to the 
exact potentials. 

As an example, we show in figure 6 both the original exact potential (broken lines), 
obtained from (]2]h and the approximate potential (solid lines) occurring after the Pekeris 
approximation, obtained from (ITT]) , for three different values of the seniority quantum 
number r, in the case of 118 Xe. In addition, the relevant energy levels, obtained from 
eq. (140 h . are shown. The following comments apply. 

1) Near their bottoms, the exact and the approximate potentials nearly coincide. 

2) The approximation breaks down, since for all values of r the relevant energy level 
lies higher than the relevant barrier on the rhs of the exact potential. We still get fits 
of good quality, which are obtained using the approximate potential, but this potential 
is not similar to the exact potential any more, preventing us to draw any conclusions 
regarding the exact potential. 

We therefore conclude that although we could derive analytical expressions for 
fitting the spectra, the parameter values coming out from the fits correspond to “well 
sizes” and diffuseness for which the approximate potential is not quite similar to the 
exact potential, i.e. the Pekeris approximation breaks down. 

6 . Conclusions 

The main results are summarized here. 

1) Approximate solutions in closed form are obtained for the 5-dimensional Bohr 
Hamiltonian with the Woods-Saxon potential, using the Pekeris approximation and the 
exact solutions of an extended Woods-Saxon potential in one dimension, featuring a 
dip near its surface. 

2 ) Applying the results to several 7 -unstable and prolate deformed nuclei, we find 
that the WS potential can describe the ground state bands and the 71 bands equally well 
as other potentials (Davidson, Kratzer, Morse), if the “well size” and the diffuseness are 
large enough (at least 1.9 and 0.44 respectively), but it fails to describe the bands, 
apparently because of its lack of a hard core. Several (forty-four) examples of deformed 
nuclei satisfying this condition have been found, but on the other hand all 7 -unstable 
nuclei considered violate this condition. 

3) The form of the potentials coming out from the fits exhibits a very large dip 
near the surface. In other words, the Bohr equation forces the parameters of the WS 
potential to obtain values producing a very large dip near its surface, so that its overall 
shape around its minimum largely resembles the shape around the minimum of the 
Davidson, or the Kratzer, or the Morse potential. 

4) The present results suggest that potentials used in the Bohr Hamiltonian can 
provide satisfactory results for nuclear spectra if they possess two features, a hard core 
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and a deep oscillator-like minimum. They also suggest that the lack of a hard core does 
not decisively affect the description of the ground state and 71 bands, but destroys the 
ability of the potential to describe /?i bands. 

Concerning the position of the quasi-/3i bands, which are not reproduced by the 
Woods-Saxon potential, the following comments apply. 

1) The quasi-/3i bandheads move to higher energies (normalized to the energy of 
the first excited state) both in 7 -unstable and in deformed nuclei, if the left wall of the 
infinite well potential used in the E(5) and X(5) critical point symmetries is gradually 
moved to the right of the origin of the /3-axis, approaching the right wall (82, [83] . 

2) The interlevel spacings within the /3i- band, which are known to be overestimated 
in the framework of the X(5) critical point symmetry, get fixed by allowing the right 
wall of the infinite well potential to be sloped to the right [84]. It should be mentioned 
here that the identification of a symmetry underlying the X(5) special solution of the 
Bohr Hamiltonian remains an open problem. 

Analytical wave functions for the Bohr Hamiltonian with the WS potential can be 
readily obtained by exploiting the similarity of this Hamiltonian, after using the Pekeris 
approximation, to the exactly soluble extended WS potential with a dip near its surface. 
The calculation of B(E2) transition rates becomes then a straightforward task, to be 
addressed in further work. 

Finally, the solution of the Bohr Hamiltonian with a Woods-Saxon potential 
within the framework of the algebraic collective model nn m mi no would offer 
the opportunity of comparison of the present approximate results to exact numerical 
solutions. 
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Table 1 . Quantum numbers appearing in the 0(5 )dS 0(3) reduction [74], occurring 
from (03 and ([48]) . 


r 

Z'A 

A 

L 

0 

0 

0 

0 

1 

0 

1 

2 

2 

0 

2 

4,2 

3 

0 

3 

6,4,3 

3 

1 

0 

0 

4 

0 

4 

8,6,5,4 

4 

1 

1 

2 

5 

0 

5 

10,8,7,6,5 

5 

1 

2 

4,2 

6 

0 

6 

12,10,9,8,7,6 

6 

1 

3 

6,4,3 

6 

2 

0 

0 
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Table 2. Comparison of theoretical predictions of the Bohr Hamiltonian with 
the Woods-Saxon potential for 46 axially symmetric prolate deformed nuclei to 
experimental data m of rare earth and actinides with -R 4/2 > 2.9 and known 0^ 
and 2+ states. The angular momenta of the highest levels of the ground state and 71 
bands included in the rms fit are labelled by L g and L 7 respectively, while N indicates 
the total number of levels involved in the fit and er is the quality measure of (1661) . All 
energies are normalized to the energy of the first excited state, E{2f). For each band, 
the i? 4/ 2 = E(4f)/E(2f) ratio (labelled by 4/2), the normalized bandhead of the 71 
band (labelled as 2 T / 2 ), and the normalized last members of the ground state and 71 
bands included in the fit (labelled by L g /2 and L 1 /2 respectively), are reported. 0o 
has been obtained from m- The theoretical predictions are obtained from (1641) . See 


subsection 0 for further discussion. 


nucleus 

Po 

A 

C 

a 0 

00 

Vo 

10 3 a o 

L g 

L ~t 

N 

0 

4/2 

exp 

4/2 

th 

Lg/ 2 
exp 

Lg/2 

th 

2 7 /2 

exp 

2 7 /2 

th 

Li/ 2 
exp 

L 7 /2 

th 

150 Nd 

4.1 

0.15 

7.6 

0.62 

0.285 

207 

43 

14 

4 

9 

0.24 

2.93 

3.02 

20.6 

20.7 

8.2 

8.6 

10.4 

9.9 

152 Sm 

4.8 

0.20 

8.6 

0.96 

0.306 

245 

61 

16 

9 

15 

0.60 

3.01 

3.10 

27.6 

28.0 

8.9 

10.2 

19.5 

18.4 

154 Sm 

4.2 

0.22 

14.0 

0.92 

0.341 

152 

75 

16 

7 

13 

0.53 

3.25 

3.26 

36.2 

36.5 

17.6 

18.6 

26.3 

25.0 

154 Gd 

3.4 

0.24 

6.9 

0.82 

0.312 

119 

75 

26 

7 

18 

0.35 

3.02 

3.14 

57.3 

57.4 

8.1 

8.8 

14.7 

14.2 

156 Gd 

2.2 

0.24 

10.8 

0.53 

0.338 

42 

81 

26 

16 

27 

0.68 

3.24 

3.25 

74.0 

74.1 

13.0 

14.5 

44.9 

44.0 

158 Gd 

2.8 

0.18 

10.6 

0.50 

0.348 

65 

63 

12 

6 

10 

0.08 

3.29 

3.28 

23.5 

23.5 

14.9 

15.1 

20.4 

20.3 

156 Dy 

3.5 

0.24 

6.0 

0.84 

0.293 

143 

70 

28 

13 

25 

0.49 

2.93 

3.08 

57.9 

57.9 

6.5 

7.5 

23.8 

22.9 

158 Dy 

2.4 

0.24 

7.7 

0.58 

0.326 

54 

78 

28 

8 

20 

0.46 

3.21 

3.21 

75.4 

74.8 

9.6 

10.3 

19.1 

18.2 

"Dy 

3.0 

0.23 

9.3 

0.69 

0.339 

78 

78 

24 

23 

33 

0.88 

3.27 

3.23 

65.1 

66.5 

11.1 

12.5 

68.2 

69.4 

160 Er 

3.0 

0.24 

5.6 

0.72 

0.304 

97 

73 

26 

5 

16 

0.52 

3.10 

3.11 

55.9 

55.4 

6.8 

7.3 

10.5 

10.0 

162 Er 

1.9 

0.24 

7.4 

0.46 

0.322 

35 

77 

20 

12 

20 

0.89 

3.23 

3.22 

43.7 

45.4 

00 

bo 

10.1 

28.5 

27.0 

164 Er 

2.1 

0.21 

6.9 
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0.333 

40 

70 
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62.8 

9.4 

9.8 

41.6 

43.1 
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83 
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4 
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164 Yb 

2.9 

0.23 

5.6 

0.67 

0.290 

100 

67 

18 

5 

12 

0.30 

3.13 

3.14 

35.6 

35.6 
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10.3 

166 Yb 
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13 
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62.3 

63.9 

9.1 
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71 

34 

7 
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64 

72 
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14 
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44 

63 
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5 

10 
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18.6 

18.8 

22.6 

22.4 

174 Yb 
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64 

65 

18 

5 

12 
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25.0 

176 Yb 

1.8 

0.20 

10.7 

0.36 

0.305 

35 

61 

18 

5 

12 

0.29 

3.31 

3.29 

48.4 

49.0 

15.4 

15.3 

19.0 

18.9 

166 Hf 

2.2 

0.25 

4.1 

0.55 

0.250 

77 

63 

22 

3 

12 

0.24 

2.97 

3.08 

36.9 

36.5 

5.1 

5.5 

6.3 

6.2 

168 Hf 

2.5 

0.24 

5.7 

0.60 

0.275 

83 

66 

22 

4 

13 

0.28 

3.11 

3.16 

46.5 

46.9 

7.1 

7.7 

9.8 

9.4 

170 Hf 

3.9 

0.22 

7.9 

0.86 

0.301 

168 

66 

34 

4 

19 

0.48 

3.19 

3.16 

105.7 

105.3 

9.5 

10.1 

12.2 

11.7 

172 Hf 

3.4 

0.23 

9.0 

0.78 

0.276 

152 

63 

38 

6 

23 

0.69 

3.25 

3.21 

132.8 

133.5 

11.3 

11.9 

17.0 

16.3 

174 Hf 

3.1 

0.23 

10.2 

0.71 

0.286 

117 

66 

22 

4 

13 

0.43 

3.27 

3.24 

58.2 

58.9 

13.5 

13.7 

15.9 

15.5 

176 Hf 

2.9 

0.22 

11.4 

0.64 

0.295 

96 

65 

18 

6 

13 

0.34 

3.28 

3.27 

45.4 

45.7 

15.2 

15.7 

21.1 

20.6 

178 Hf 

2.2 

0.22 

9.2 

0.48 

0.280 

62 

62 

18 

6 

13 

0.41 

3.29 

3.26 

44.2 

45.0 

12.6 

12.9 

18.1 

17.8 

176 W 

2.4 

0.24 

7.6 

0.58 




22 

5 

14 

0.31 

3.22 

3.21 

51.8 

51.8 

9.6 

10.2 

14.0 

13.2 

178 W 

2.4 

0.23 

7.6 

0.55 




14 

2 

7 

0.14 

3.24 

3.23 

27.0 

27.2 

10.5 

10.4 

10.5 

10.4 

180 w 

2.7 

0.24 

8.7 

0.65 

0.254 

113 

61 

24 

7 

17 

0.67 

3.26 

3.22 

60.0 

60.4 

10.8 

11.6 

18.7 

17.6 

176 Os 

4.0 

0.23 

6.0 

0.92 




24 

5 

15 

0.37 

2.93 

3.04 

45.5 

45.5 

6.4 

7.3 

10.4 

9.6 

178 Os 

3.8 

0.23 

6.0 

0.87 




16 

5 

11 

0.40 

3.02 

3.06 

26.1 

26.0 

6.6 

7.4 

10.8 

9.9 

180 Os 

2.5 

0.25 

5.9 

0.63 

0.226 

122 

57 

14 

7 

12 

0.59 

3.09 

3.14 

21.8 

22.0 

6.6 

7.7 

14.2 

13.0 

182 Os 

3.7 

0.24 

6.7 

0.89 

0.234 

250 

56 

26 

7 

18 

0.84 

3.15 

3.10 

54.0 

53.3 

7.0 

8.3 

14.6 

13.4 

184 Os 

2.2 

0.24 

6.2 

0.53 

0.213 

107 

51 

22 

6 

15 

1.05 

3.20 

3.19 

47.9 

49.4 

7.9 

8.4 

13.5 

12.8 

228 Ra 

3.1 

0.24 

10.1 

0.74 

0.217 

204 

52 

22 

3 

12 

0.18 

3.21 

3.23 

53.6 

54.0 

13.3 

13.3 

14.1 

14.0 

22 S 

3.5 

0.24 

13.1 

0.84 

0.230 

231 

55 

18 

5 

12 

0.15 

3.24 

3.25 

41.7 

41.6 

16.8 

17.1 

20.3 

20.0 

23 U 

2.8 

0.22 

10.6 

0.62 

0.244 

132 

54 

18 

4 

11 

0.12 

3.27 

3.26 

45.1 

45.3 

14.7 

14.7 

16.6 

16.6 

232 

2.4 

0.23 

12.1 

0.55 

0.261 

85 

60 

30 

12 

25 

0.54 

3.28 

3.27 

104.6 

105.9 

15.9 

16.6 

36.5 

35.6 

232 u 

3.0 

0.22 

13.3 

0.66 

0.264 

129 

58 

20 

4 

12 

0.20 

3.29 

3.28 

55.9 

56.3 

18.2 

18.4 

20.4 

20.3 

234 u 

2.7 

0.23 

15.9 

0.62 

0.272 

99 

63 

28 

7 

19 

0.36 

3.30 

3.29 

98.8 

99.6 

21.3 

21.8 

29.0 

28.5 

236 

3.3 

0.22 

15.5 

0.73 

0.282 

137 

62 

26 

5 

16 

0.63 

3.30 

3.29 

89.3 

90.6 

21.2 

21.4 

24.9 

24.7 

238 

2.9 

0.23 

18.3 

0.67 

0.286 

103 

66 

30 

27 

40 

0.85 

3.30 

3.29 

112.1 

114.9 

23.6 

25.0 

112.7 

113.5 

24 »Pu 

2.2 

0.22 

19.0 

0.48 

0.289 

58 

64 

26 

4 

15 

0.21 

3.31 

3.30 

95.5 

96.0 

26.6 

26.6 

28.8 

28.6 

242 Pu 

2.3 

0.22 

17.7 

0.51 

0.292 

62 

64 

26 

2 

13 

0.51 

3.31 

3.30 

93.7 

94.8 

24.7 

24.8 

24.7 

24.8 

248 Cm 

2.7 

0.22 

17.2 

0.59 

0.297 

83 

65 

28 

2 

14 

0.87 

3.31 

3.30 

10.5 

10.7 

24.2 

24.0 

24.2 

24.0 
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Table 3. Spectra of some rare earths and actinides, obtained with the parameters 
shown in Table [2j 


L 

156 Gd 

exp 

156 Gd 

th 

164 Er 

exp 

164 Er 

th 

17°Y b 

exp 

17°Y b 

th 

232 Th 

exp 

232 Th 

th 

gsb 









4 

3.24 

3.25 

3.28 

3.25 

3.29 

3.26 

3.28 

3.27 

6 

6.57 

6.59 

6.72 

6.60 

6.80 

6.66 

6.75 

6.69 

8 

10.85 

10.86 

11.21 

10.92 

11.43 

11.07 

11.28 

11.13 

10 

15.92 

15.88 

16.61 

16.11 

17.06 

16.37 

16.75 

16.47 

12 

21.63 

21.56 

22.79 

22.09 

23.54 

22.49 

23.03 

22.60 

14 

27.83 

27.78 

29.58 

28.83 

30.63 

29.36 

30.04 

29.45 

16 

34.39 

34.50 

37.33 

36.29 

37.92 

36.94 

37.65 

36.97 

18 

41.29 

41.66 

45.10 

44.44 

45.18 

45.21 

45.84 

45.11 

20 

48.62 

49.23 

53.28 

53.27 

52.66 

54.13 

54.52 

53.85 

22 

56.49 

57.18 

61.85 

62.77 



63.69 

63.16 

24 

64.95 

65.49 





73.32 

73.03 

26 

73.99 

74.15 





83.38 

83.45 

28 







93.82 

94.39 

30 







104.56 

105.86 

7i 









2 

12.97 

14.49 

9.41 

9.80 

13.60 

14.02 

15.91 

16.59 

3 

14.03 

15.27 

10.36 

10.65 

14.54 

14.85 

16.80 

17.41 

4 

15.24 

16.29 

11.58 

11.77 

15.78 

15.95 

18.03 

18.49 

5 

16.94 

17.55 

13.10 

13.15 

17.33 

17.31 

19.45 

19.82 

6 

18.47 

19.04 

14.87 

14.78 

19.01 

18.92 

21.27 

21.41 

7 

20.79 

20.74 

16.91 

16.65 

21.13 

20.78 

23.21 

23.23 

8 

22.61 

22.65 

19.09 

18.76 

23.19 

22.87 

25.50 

25.29 

9 

25.29 

24.74 

21.64 

21.08 

25.76 

25.18 

27.75 

27.57 

10 

27.45 

27.02 

23.90 

23.63 

28.16 

27.72 

30.62 

30.06 

11 

30.20 

29.47 

27.13 

26.38 

30.90 

30.47 

33.22 

32.75 

12 

32.85 

32.08 

29.91 

29.33 

33.55 

33.42 

36.48 

35.64 

13 

35.69 

34.84 

33.13 

32.48 

36.40 

36.57 



14 

38.64 

37.74 

35.75 

35.83 

39.25 

39.91 



15 

41.76 

40.78 

38.51 

39.36 





16 

44.90 

43.95 

41.59 

43.07 
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Table 4. Same as Table O for 17 additional nuclei. See subsection 15.11 for further 
discussion. 


nucleus 

00 

A 

C 

a 0 

00 

Vo 

10 3 a 0 

L 9 

L-y 

N 

£7 

4/2 

exp 

4/2 

th 

L g /2 

exp 

Lg/2 

th 

2 7 /2 

exp 

2 7 /2 

th 

£ 7 / 2 
exp 

1 . 7/2 

th 

160 Gd 

1.8 

0.18 

9.0 

0.32 

0.353 

26 

64 

16 

8 

14 

0.10 

3.30 

3.29 

40.0 

40.2 

13.1 

13.1 

22.8 

22.8 

162 Gd 

1.5 

0.18 

8.2 

0.27 




14 

4 

9 

0.05 

3.29 

3.29 

31.4 

31.5 

12.0 

12.0 

14.1 

14.1 

162 Dy 

1.9 

0.19 

7.7 

0.36 

0.343 

31 

65 

18 

14 

21 

0.23 

3.29 

3.28 

47.6 

47.9 

11.0 

11.1 

39.4 

39.2 

164 Dy 

1.2 

0.19 

7.1 

0.23 

0.348 

12 

66 

20 

10 

18 

0.26 

3.30 

3.28 

57.4 

58.1 

10.4 

10.4 

25.3 

25.3 

166 Dy 

1.2 

0.15 

7.5 

0.18 




6 

5 

6 

0.02 

3.31 

3.30 

6.9 

6.8 

11.2 

11.2 

14.9 

14.9 

166 Er 

1.2 

0.21 

6.9 

0.25 

0.342 

12 

72 

16 

14 

20 

0.30 

3.29 

3.26 

36.8 

37.1 

9.8 

9.9 

35.7 

36.6 

!68 B r 

1.2 

0.16 

6.9 

0.19 

0.338 

13 

54 

18 

8 

15 

0.20 

3.31 

3.30 

50.0 

50.6 

10.3 

10.3 

20.4 

20.3 

170 Br 

0.7 

0.18 

8.9 

0.13 

0.336 

4 

61 

26 

19 

30 

0.85 

3.31 

3.30 

95.8 

98.1 

11.9 

13.0 

66.2 

65.7 

178 Yb 

1.3 

0.17 

9.9 

0.22 




6 

2 

3 

0.01 

3.31 

3.30 

6.9 

6.9 

14.5 

14.5 

14.5 

14.5 

180 H f 

1.0 

0.16 

8.7 

0.16 

0.274 

13 

44 

12 

5 

9 

0.05 

3.31 

3.30 

24.3 

24.4 

12.9 

12.9 

16.7 

16.6 

182 w 

1.4 

0.20 

8.6 

0.28 

0.251 

31 

50 

18 

6 

13 

0.26 

3.29 

3.28 

47.4 

48.0 

12.2 

12.4 

17.7 

17.6 

I8 4 w 

0.9 

0.18 

5.4 

0.16 

0.236 

15 

43 

10 

6 

9 

0.06 

3.27 

3.28 

16.7 

16.8 

8.1 

8.0 

13.3 

13.4 

186 w 

1.0 

0.19 

4.1 

0.19 

0.226 

20 

43 

14 

6 

11 

0.09 

3.23 

3.25 

29.1 

29.2 

6.0 

6.1 

11.4 

11.4 

186 Os 

1.7 

0.22 

4.2 

0.37 

0.200 

72 

44 

14 

13 

18 

0.20 

3.17 

3.19 

25.9 

26.1 

5.6 

6.1 

26.5 

26.8 

188 Os 

1.4 

0.23 

3.0 

0.32 

0.186 

57 

43 

12 

7 

11 

0.18 

3.08 

3.14 

18.4 

18.7 

4.1 

4.4 

10.9 

10.6 

190 Os 

1.4 

0.24 

2.2 

0.34 

0.178 

62 

43 

10 

6 

9 

0.20 

2.93 

3.06 

12.6 

12.7 

3.0 

3.3 

7.9 

7.6 

288 Pi] 

1.7 

0.21 

16.5 

0.36 

0.286 

35 

60 

26 

4 

15 

0.34 

3.31 

3.30 

96.8 

97.6 

23.3 

23.4 

25.5 

25.5 


Table 5. Comparison of theoretical predictions of the Bohr Hamiltonian with the 
Woods-Saxon potential to experimental data EH for 9 7 -unstable Xe isotopes with 
i? 4/ 2 < 2.6 and known Of and 2+ states. The angular momenta of the highest levels of 
the ground state and quasi -71 bands included in the rms fit are labelled by L g and L 7 
respectively, while N indicates the total number of levels involved in the fit and a is the 
quality measure of (1661) . All energies are normalized to the energy of the first excited 
state, E(2f). For each band, the i ? 4 / 2 = E(4:f)/E(2f) ratio (labelled by 4/2), the 
normalized bandhead of the quasi -71 band (labelled as 2 7 / 2 ), and the normalized last 
members of the ground state and quasi -71 bands included in the fit (labelled by L g /2 
and L 7 /2 respectively), are reported. 0q has been obtained from [71]. The theoretical 
predictions are obtained from (l44l) . See subsection 15.21 for further discussion. 


nucleus /3 q A 10 3 ao Vo 10 3 ao L g LN u 4/2 4/2 L g /2 L g /2 2-y/2 2-y/2 L-y/2 L.y/2 

exp th exp th exp th exp th 


118 Xe 

0.13 

0.18 

23 

0.265 

0.2 

48 

14 

8 

13 

120 Xe 

0.21 

0.17 

36 

0.291 

0.5 

50 

14 

9 

14 

122 Xe 

0.52 

0.14 

71 

0.259 

4.0 

35 

16 

9 

15 

124 Xe 

0.57 

0.16 

88 

0.212 

7.2 

33 

12 

8 

12 

126 Xe 

0.65 

0.18 

115 

0.188 

11.9 

33 

12 

9 

13 

128 Xe 

0.67 

0.21 

143 

0.184 

13.4 

39 

10 

7 

10 

130 Xe 

0.92 

0.22 

203 

0.169 

29.4 

37 

14 

5 

10 

132 Xe 

0.90 

0.24 

218 

0.141 

41.2 

34 

6 

2 

3 

134 Xe 

1.10 

0.27 

302 

0.119 

85.6 

33 

8 

5 

7 


0.26 

2.40 

2.34 

12.9 

13.1 

2.8 

2.3 

7.8 

8.0 

0.47 

2.47 

2.37 

13.8 

13.6 

2.7 

2.4 

9.8 

10.8 

0.65 

2.50 

2.32 

17.5 

17.6 

2.5 

2.3 

9.7 

11.1 

0.51 

2.48 

2.32 

11.0 

10.8 

2.4 

2.3 

8.2 

8.2 

0.63 

2.42 

2.28 

11.0 

10.1 

2.3 

2.3 

9.1 

10.1 

0.47 

2.33 

2.20 

7.6 

6.9 

2.2 

2.2 

6.2 

6.9 

0.30 

2.25 

2.11 

9.5 

9.5 

2.1 

2.1 

4.1 

4.7 

0.34 

2.16 

2.05 

3.2 

3.2 

1.9 

2.1 

1.9 

2.1 

0.29 

2.04 

1.82 

3.5 

3.1 

1.9 

1.8 

2.7 

3.1 








